F/6  9/5 


AO-A047  411  naval  RESEARCH  LAB  MASHIN6T0N  0 C 

VOLTAGE-CONTROLLED  OSCILLATOR  WITH  TIME-DELAY  FEEDBACK. (U) 
SEP  77  J B RAO»  W M WATERS 

UNCLASSIFIED  NRL-8138  SBIE-AD-E000031 


NL 


I 

' I 


1 i 


AD  A 0 47^ 


September  8,  1977 


Apfiroved  for  puNic  icleuc;  dinributioa  imlimiud. 


SKCu'tTv  Cl  ASSlFlCAT'ON  Of  TmiS  PAGE  ^'»l‘h«n  Data  f.nfr^dj 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

' REPORT  NLINSER 2 GOVT  ACCESSION  NO. 

f NRL~BP|IH138  / 

RECIPIENT'S  CATAlOCs  NUMBER 

4 Title  ^And  Subf/fl*)  . 

^4.  Type  AP  RfonRT  « PERIOD  COVERED 

j J/OLTAGE-CONTROLLED  (g^SCILLATOR  WITH  ■ 

Interim  /tep^rtpn  a continuing 

NRL  Problem 

TIME-DELAY  FEEDBACK  § - 

- t?  JA 

6 PERFORMING  ORG.  REPORT  NUMBER 

B.  contract  or  grant  NUMBCRf*; 

y J.  B.  L^Rao  W.  M/waters  j 

7^/x  /ei  \ 

9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Research  Laboratory 

Washington,  D.  C.  20375 

f]} 

10  program  element  PROJECT.  TASK 
AREA  & WORK  UNIT  NUMBERS 

NRL  Problen^RQ2-92^01 

Projecl^  Z^12J5f  001  | 

u.  controlling  OFFICE  NAME  And  ADDRESS  ' H 

Naval  Research  Laboratory  ' //  ) 

-u — Mooaj 

Washington,  D.  C.  20375  /J' 

<3  NUMBER  OF  PAGES 

14 

14  monitoring  agency  NAME  A AODRESSCIf  dtfforont  from  Conirolling  OIttco) 

15.  SECURITY  CLASS,  (ot  roport. 

Unclassified 

ISa  OECL  AbSiFiCATlON  DOWNCRAO'NG 
SCHEDULE 

16  Distribution  statement  (cI  this  Roport) 

Approved  for  public  release;  distribution  unlimited. 

/V  - - J //J/-  rps- . ' ' -■  / 

>7  distribution  statement  fo/  fha  abaffacf  •nl0red  In  Block  20.  H dm9ront  trom  Roport) 


I®  suPPL ementarv  notes 


t9  KEY  WORDS  fConilnum  on  ravaraa  oido  if  naraaaao  Idonflly  6^  block  numbot) 

Circuits  Radar 

Coherence 

Microwaves 

Oscillator 


tract  ^Conrlnua  on  ravaraa  alda  If  nacaaaary  and  idaniUy  by  block  numbor) 
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VOLTAGE-CONTROLLED  OSCILLATOR  WITH 
TIME-DELAY  FEEDBACK 


INTRODUCTION 

A voltage-controlled  oscillator  (VCO)  with  time-delay  feedback  has  been  designed  as 
a frequency-agile,  noise-degenerated  radar  frequency  source.  This  frequency  source  uses  a 
solid  state  VCO  with  a noise-degenerated  negative  feedback  loop  with  an  interferometer  as 
a frequency  discriminator.  Because  of  the  periodicity  of  responses  of  the  interferometer, 
the  oscillator  operation  is  stable  at  many  frequencies  across  the  band  of  interest.  For  this 
reason,  both  low  noise  from  the  degeneration  and  frequency  agility  from  the  periodic 
response  are  obtained  in  one  source.  This  report  analyzes  the  steady  state  and  transient 
responses  of  the  frequency  source.  The  steady  state  solution  will  provide  information 
necessary  in  choosing  the  loop  parameters  that  set  the  noise  degeneration  level.  The  trans- 
ient response  yields  information  on  loop  response  time,  and  thus  the  switching  time  needed 
for  frequency  agility. 

GENERAL  DESCRIPTION 

Figure  1 is  a block  diagram  of  a voltage-controlled  oscillator  (VCO)  with  a time-delay 
feedback  loop  containing  an  interferometer  and  a video  amplifier.  The  interferometer’s 
output  is  proportional  to  the  phase  difference  between  the  direct  and  delayed  outputs 
from  the  VCO.  If  the  VCO  output  is  given  by  cos(cuO.  the  interferometer  output  after 
the  high-frequency  component  is  filtered  out,  is  given  by 

e(f)  = sincof. 


Fig.  1 — Voltage-controlled  oscillator  with  time-delay  feedback 
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Fig.  2 — Interferometer  response  with 
additional  phase  shift 


Therefore,  interferometer  output  voltage  vs  frequency  is  a sine  wave,  so  responses  repeat, 
spaced  by  a frequency  equal  to  the  inverse  of  the  delay  time,  as  shown  in  Fig.  2.  For  a 
delay  time  of  25  ns,  the  periodicity  equals  40  MHz.  The  whole  response  curve  can  be 
shifted  by  adding  additional  phase  shift  in  the  direct  path  [IJ.  A 90-degree  phase  shift 
will  shift  the  whole  response  pattern  one-fourth  of  its  period  (10  MHz).  In  this  manner, 
the  interferometer  can  operate  with  a ±20-MHz  pull-in  range  and  still  produce  a response 
every  10  MHz  across  the  band.  Because  of  this  periodicity  of  responses,  stable  operation 
occurs  at  many  frequencies  across  the  band  of  interest,  providing  frequency  agility. 


ANALYSIS 

The  system  shown  in  Fig.  1 was  analyzed,  with  the  input  voltage  as  a step  function 
of  an  amplitude  e,  to  study  the  effectiveness  of  the  time  delay  feedback.  When  input 
voltage  V is  zero,  the  output  angular  frequency  of  the  VCO  is  assumed  to  be  and  is 
related  to  time  delay  r as  = 2irN/T,  where  N is  an  integer.  When  this  relation  holds, 
becomes  one  of  the  stable  frequencies,  and  the  output  voltage  of  the  interferometer 
will  be  zero.  Now,  if  input  voltage  v is  applied,  the  output  frequency  of  the  VCO  will 
change.  That  change  in  frequency  is  given  by 


Aw  = -K2  e^(t) 


(1) 
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where  6^(1)  is  the  voltage  at  the  VCO  input  and  K2  is  the  modulation  sensitivity  in  radians 
per  volt. 

Let  0(t)  represent  the  phase  shift  due  to  this  frequency  change  during  a time  t.  The 
quantity  <t>(t)  is  given  by 


m = - K2  e,(t)  dt.  (2) 

*^o 

Then,  interferometer  output  voltage  e(t),  after  the  high-frequency  component  is  filtered 
out,  is 

e(t)  = Kj  sin[0(t)  - 0((  - t)],  (3) 

where  Kj  is  the  phase-detector  sensitivity  constant  in  volts  per  radian. 

For  [(/)(f)  - 0(t  - t)]  < 1,  which  is  satisfied  near  the  stable  frequency  of  operation, 
one  can  approximate: 


e{t)  ^ Klim  - <P(t  - T)]. 

(4) 

Taking  Laplace  transforms  of  Eqs.  (1),  (2),  and  (4)  yields 

An(s)  = -K2  Ec(s) 

(5) 

<I>(s)  = -K2  Emis 

(6) 

E{s)  = Ki[‘I>(s)  -4>(s)e-^»]. 

(7) 

In  addition.  From  Fig.  1 we  have 

E,(s)  =-  + £(8)  G(s), 
s 

(8) 

where  G(s)  is  the  transfer  function  of  the  video  amplifier. 

Using  Eqs.  (5)  to  (8),  we  can  show  that 

-Kov 

Af2(s) 

KlX2G(s)[l-e-’'«]  +s 

(9) 

The  amplifier  in  the  feedback  loop  is  assumed  to  be  a combination  video  amplifier 
and  low-pass  filter  with  combined  transfer  function 

a G„ 

G{8) 

8 + a 

(10) 

where  is  the  low-frequency  gain  of  the  amplifier  and  a is  the  3-dB  cutoff  point. 
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Substituting  Eq.  (10)  in  Eq.  (9)  gives  the  normalized  response 


F(s) 


s + a 


An(s)  _ 

-K2V  g2  + sa  + i4(l  - e"’"*) 


(11) 


where  A = /fj  K2  do®- 


Steady  State  Solution 

A steady  state  solution  can  be  obtained  easily  from  Eq.  (11),  as 


Aco(<») 

- ^ = Urn  AO  = Lim  sF(s)  = 

-K2I’  / -*  00  s ^ o 


1 

1^2  Go  1 


(12) 


This  steady  state  solution  can  be  used  [2]  to  find  the  amount  of  VCO  noise  degenera- 
tion possible  with  a time-delay  feedback  system.  Let  represent  the  PM  noise  of 

VCO  without  feedback.  The  corresponding  equivalent  noise  voltage  at  the  VCO  input  is 
27rAA,s(./(-F2).  Let  A/  be  the  VCO  FM  noise  with  feedback.  Then,  substituting 
i'  = 2iTAf„g^f(-K2)  and  A<jj(°°)  = 27r A f in  Eq.  (12),  we  obtain 


Af  = 


K,  Ko  r + I 


(13) 


From  Eq.  (13)  it  may  be  noted  that  the  VCO  noise  is  degenerated  by  a factor  of 
(Ki  K2  Go  T + 1).  Therefore,  to  reduce  VCO  noise,  it  is  necessary  that 

(Kl  K2  Go  r + l)>l.  (14) 


Stability  Condition 

The  condition  for  stability  of  the  feedback  system  is  discussed  in  Appendix  A and  is 
given  by 


Kj  K2  Go«t2  <7r2/2. 


(15) 


which  is  identical  to  Eq.  (All). 

For  the  system  to  be  stable  and  produce  noise  degeneration,  K2,  Gq,  a,  and  t 
should  be  chosen  such  that  the  conditions  in  both  Eqs.  (14)  and  (15)  are  satisfied. 


Transient  Response 

Formally,  the  inverse  Laplace  transform  of  Bq.  (11)  gives  the  transient  response. 
However,  the  inverse  transform  is  not  readily  available.  The  following  procedure  [3]  is 
used  to  obtain  the  transient  response.  Repeating  Eq.  (11),  we  have 
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S + tt 

s'^  + set  + A - 


Let  z = + sa  + ^ . Then 


If  the  expansion  of  the  type 


is  used,  Eq.  (17)  becomes 


F(s) 


s + a 1 


1 - X 


^ 1 ^ -TS 

1 — c '* 
z 


1 + or  + 


(16) 


(17) 


F(s)  = 1 1 + — e-^*  + — c-2«  + ...I . 


(18) 


Taking  the  inverse  Laplace  transform  of  Eq.  (18)  yields 

f(t)  = /j(()  + Au{t  - T)f2(t  - t)  + A^u(t  - 2T)f^{t  - 2t)  + ... 

where 

fn(t)  = I'-*  r --  ^ 1 

L(s2  + SO  + .4)"J 

and  u{t)  is  a unit  step  function. 


(19) 


(20) 


The  transient  response  given  in  Eq.  (19)  appears  to  contain  an  infinite  number  of 
terms.  However,  for  any  given  finite  time  f,  only  a finite  number  of  terms  are  nonzero. 
Also,  each  individual  term  has  a physical  significance  related  to  the  time-delay  feedback. 

In  the  time  range  0 < t < r,  only  the  first  term  is  nonzero  and  the  rest  of  the  terms  are 
zero,  so  that  the  first  term  gives  the  system  response  before  the  VCX)  output  is  applied 
to  the  phase  detector  through  the  delayed  path.  In  the  range  t < ( < 2t,  the  signal 
through  the  delayed  path  is  applied  to  the  phase  detector,  and  the  transient  response  is 
represented  by  the  first  two  terms,  since  the  remaining  terms  are  still  zero.  In  general, 
for  range  (Af  - 1)t  < t < Aft,  the  first  Af  terms  give  the  transient  response  because  the 
higher  order  terms  are  still  zero. 

The  transient  response  is  not  yet  complete,  because  the  inverse  transform  shown  in 
Eq.  (20)  is  not  readily  available.  A closed  form  solution  for  f„(t)  is  obtained  by  using  the 
following  procedure.  First,  consider  the  inverse  Laplace  transform. 


S„(t)  = L-l 


s2  + sa  + A)” 


This  can  be  written  as 
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SJt)  = L-1 


(s  + a)"  (s  + bY‘ 


(21) 


where  o = — - joji . ^ “ IT  . ^nd  cjj  = yA  - ——.  It  will  be  shown  later  that  a 2/4 

iL  £t  f *^1 

Therefore,  toj  is  real. 


by 


The  inverse  Laplace  transform  shown  in  Eq.  (21)  is  readily  available  [4].  It  is  given 


_ v/^ 


(n  - 1)!  \a  - b 


ii-'A 


{a  + b)t 


‘n-'A 


a + b\ 


where  /„_i^  is  a modified  Bessel  function  of  the  first  kind.  Substituting  for  a and  b from 
Eq.  (21),  we  obtain 


S„(0  = — 


{n  - 1)!  V2co 


a/ 


c 2 (wjO 


(22) 


where  is  the  Bessel  function  of  the  first  kind. 


From  Eqs.  (20)  and  (21)  it  may  be  noted  that  which  can  be  obtained  if  s„(f) 
is  known,  is  given  by 


/'n(0  = s„(0  - s„(+0)  + Oi  s„(0 
where  s„(f)  is  the  first  derivative  of  .Sn(f)- 


(23) 


By  noting  s„(+0)  = 0 and  substituting  s„(t)  from  Eqs.  (22)  and  (23)  and  simplifying, 
we  can  show  that 


v/tT  / t 
” (n  - 1)!  (^1 


n-'A 


at 


^‘*^1  *^0-3/2 0 +"2 ‘^n-l/2(‘^l^) 

In  terms  of  spherical  Bessel  functions  j^,  one  can  show  [5]  that 


(24) 


/„(0  = 


(n 


n-l 


at  _ 


(cJiO  e 


;„.2(cjjf)  +-^  ;„-i(wiO 


(25) 


Hence,  the  transient  response  is  given  by  Eq.  (19),  with  /■„(!)  as  shown  in  Eq.  (25). 
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NUMERICAL  RESULTS 

With  Eqs.  (19)  and  (25),  the  transient  response  is  computed  for  typical  values  of  the 
system  parameters,  = 0.2  V/rad,  = 3.5  X 10®  rad/sV,  a = 6.283  X 10®  rad/s 
(corresponds  to  an  amplifier  filter  bandwidth  of  1 kHz),  and  t = 25  ns.  Figure  3 shows 
the  normalized  transient  response  with  the  above  parameters  and  for  the  amplifier  gains 


Fig.  3 — Transient  response 

of  60  dB  (G„  = 10®),  70  dB  (G„  = 3.162  X 10®),  80  dB  (G^  = 10^),  85.0823  dB  (G^  = 
1.7952  X 10'*),  and  90  dB  (G^,  = 3.162  X 10"*).  The  corresponding  values  of  At® 

( = KjX2G„ar2)  are  0.275,  0.869,  2.75,  7r2/2,  and  8.69.  For  the  system  to  be  stable. 

At®  must  be  less  than  n^/2  ( = 4.9348).  From  Fig.  3,  it  is  evident  that  for  Ar^  > 

4.9348  the  system  transient  response  is  oscillatory  with  increasing  amplitude  with  time.  This 
indicates  that  the  system  is  unstable.  For  At®  = 7r®/2,  the  response  is  oscillatory  with 
approximately  constant  amplitude,  indicating  that  the  stability  condition  given  by  Eq.  (15)  is 
highly  accurate.  For  At®  = 2.75,  the  transient  response  is  oscillatory,  but  the  oscillations 
are  damped  with  time,  suggesting  that  the  system  is  stable.  As  the  amplifier  gain  is  reduced, 
making  At®  smaller,  the  damped  oscillations  almost  vanish  for  At®  = 0.869  (which  may  be 
called  the  critically  damped  case).  For  smaller  values  of  At®,  the  transient  response  decreases 
with  time  (similar  to  an  exponentially  decreasing  curve)  and  takes  longer  to  reach  steady  state. 

These  results  indicate  that  a value  for  At®  =»  1 may  be  a best  compromise  for  faster 
loop  response  and  acceptable  damped  oscillations.  For  the  critically  damped  case,  the 
loop  response  reaches  steady  state  in  a time  span  of  about  8t.  For  a given  value  of  At®, 
the  loop  response  time  is  smaller  for  smaller  values  of  t,  as  shown  in  Fig.  4.  However, 
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Fi^.  4 — Frequency  response  wilh  T as  a parameter 


if  the  time  is  normalized  with  respect  to  r,  the  transient  responses  for  both  cases  are 
approximately  similar,  except  that  the  steady  state  respon'  is  smaller  for  smaller  r.  Since 
,4  = A'jA'gGoQ,  changing  and  A'g  will  have  the  same  effe  ' as  changing  However, 
changing  a has  approximately  the  same  effect  on  the  transi-  ..t  response  as  changing  G„, 
but  the  steady  state  response  is  independent  of  cv.  For  these  rea-  )ns,  a should  be  chosen 
as  small  as  possible,  and  KjK2C„  should  be  chosen  as  large  as  possible,  for  a given  t such 
that  At^  = 1 and  Aj  A2G„r  > 1,  which  will  satisfy  the  requirements  of  noise  degeneration 
and  system  stability  conditions  given  by  Eqs.  (14)  and  (15). 

CONCLUSIONS 

A voltage-c'ontrolled  oscillator  with  time-delay  feedback  has  been  analyzed.  By  properly 
choosing  system  piu'ameters,  low  noise  from  the  degeneration  and  frequency  agility  from  the 
periodic  response  can  be  obtained  in  one  source.  The  conditions  for  system  stability  and  the 
amount  of  noise  degeneration  are  expressed  in  terms  of  system  parameters.  Effects  of  differ- 
ent parameters  on  system  stability  and  noise  degeneration  are  discussed  in  detail.  It  is  shown 
that  the  open-loop  gain  should  be  chosen  as  large  as  possible  for  good  noise  degeneration, 
and  the  bandwidth  should  be  adjusted  to  meet  the  stability  condition. 
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Appendix  A 

STABILITY  CRITERION 


The  Nyquist  stability  criterion  is  used  to  find  the  condition  under  >xh,ch  the  feedback 
system  is  stable  and  to  identify  the  system  parameters  on  which  this  condition  de,H>n  s 
This  criterion  furnishes  a graphical  method  for  determining  the  stability  of  a sysU  m.  The 
system  is  stable  if  a plot  of  its  open-loop  transfer  function  Z(s)  for  a succession  of  values 
of  s encircling  the  entire  right  half  of  the  s-plane  in  the  clockwise  direction,  makes  a 
number  of  counterclockwise  revolutions  about  the  critical  point  (-1  + ;0)  equal  to  the 
number  of  the  poles  of  Z(s)  in  the  right  half  of  the  s-plane. 


For  the  system  in  Fig.  Al,  the  open-loop  transfer  function  is  given  by 


Z{s) 


A(l-c-^^) 
s(s  + a) 


(Al) 


From  Eq.  (Al),  note  that  there  are  no  poles  in  the  right  half  of  the  s-plane.  The  only 
pole  at  S = -a  is  in  the  left  half.  Also,  there  are  no  poles  on  the  imaginary  axis.  Since 
7(sl  becomes  zero  for  s to  apply  the  Nyquist  criterion  one  needs  a polar  plot  of 
zls!  withT=  /CO  and  for  the ’range’  of  values  - < co  < « only.  As  can  be  noted  a polar 
plot  for  the  frequency  range  -0°  to  0 is  the  mirror  im^e  about  the  ^ 

the  plot  for  the  frequency  range  0 to  «>.  Thus,  to  evaluate  the  polar  plot  of  ^ transfer 
function,  the  only  frequency  range  to  be  considered  is  from  0 to 

s = /CO  in  Eq.  (Al)  and  finding  real  part  Zr(;co)  and  imaginary  part  Z,(;(o)  of  transfer 
function  Z(/co),  it  can  be  shown  that 


Zii  (/co)  = 


a2  + 


' sincoT  „ . 2 cor 

(XT 2 sin^  — 

cor  2 , 


(A2) 


2/(/co)  = 


-A 


(X^  + co^ 


sincor  + ar- 


sin(cor/2)  sin(cor/2) 
(cor/2) 


(A3) 
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From  Eqs.  (A2)  and  (A3),  note  that  there  is  no  simple  way  of  sketching  the  polar 
plot  if  parameters  A,  a,  and  t are  not  specified.  Fortunately,  for  the  system  under  con- 
sideration the  range  of  parameters  is  such  that  a simplification  is  possible.  For  at  of  the 
order  of  10”^,  the  first  term  in  the  parentheses  on  the  right  hand  side  of  Eq.  (A2)  can  be 
neglected,  except  possibly  when  cot  =*=  2Nrr.  (In  that  case,  both  terms  will  he  dost*  to 
zero  and  Z|j(;co)  will  be  very  small.)  Similarly,  the  second  term  in  the  brackets  of  Eq.  (A3) 
can  be  neglected,  except  possibly  at  or  very  near  cor  = Ntt.  In  that  cast",  both  terms  will  be 
very  small,  which  in  turn  makes  Z/  small.  These  approximations  and  observations  are  used 
only  to  simplify  the  procedure  in  obtaining  the  polar  plot.  No  approximation  is  needed  in 
obtaining  the  exact  condition  for  stability,  as  will  be  shown  later.  K(*eping  aside  the  term 
A/{a^  + cj2),  which  is  common  for  both  Z^j  and  Z/,  it  may  be  noted  that  Z^  varies  as 
-2  sin2( cot/2)  and  Z/  varies  as  -sin(coT).  These  two  functions  are  sketched  in  Fig.  A2  as  a 
function  of  cot. 


Fit>.  A2  — Variation  of  and  except 
for  the  factor  AI(oZ  + cO*) 


In  the  range  of  parameters  of  interest,  cv  is  of  the  order  of  10“^  and  t is  of  the  order 
of  10"*^.  Therefore,  for  cot  > 7r/2,  co  > a and  the  factor  A/(a2  + co^)  = A/co^.  Using 
these  observations  and  Fig.  A2,  one  can  sketch  the  polar  plot  of  Z(;co),  which  is  shown 
in  Fig.  A3  for  the  range  0 < co  < o®.  For  the  range  -o®  < co  < 0,  the  polar  diagram  will 
be  the  mirror  image  about  the  Z^  axis,  and  it  is  not  shown.  From  the  sketch  in  Fig.  A3 


Fig.  A3  — Sketch  of  a polar  plot  for 
Z(yto),  for  A > Ap 


it  is  clear  that  if  the  real  part  Zy;  < -1  ( IZyj  I >l)  when  Z;  = 0,  the  critical  point  is 
encircled  counterclockwise  by  the  polar  plot.  Therefore,  the  system  ^ill  be  unstable.  Since 
the  magnitude  of  Z/j  depends  on  A,  for  the  system  to  be  stable  A should  be  smaller  than 
that  for  which  Zyj  = -1  and  Zj  = 0. 

Next,  it  is  possible  to  find  the  value  of  A corresponding  to  Zyj  = -1  and  Z/  = 0.  This 
will  be  done  using  Eqs.  (A2)  and  (A3),  first  without  making  any  approximations.  Equating 
Zf  = 0,  we  obtain  from  Eq.  (A3) 
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sinwT  = -2  — sin-^  — 
a;  2 


(A4) 


Equating  = -1,  we  obtain  from  Eq.  (A2) 


A / « . „ . 9 CUT)  , 

— sincoT  - 2sin^  — = -1. 

Cv2  + \ CJ  2 j 


(A5) 


Substituting  from  Eq.  (A4)  for  sintor  in  Eq.  (A5)  yields 


„ . . COT  , 

2.4  sm^  — = co'‘. 


(A6) 


From  Eqs.  (A4)  and  (A6),  we  obtain 


cos 


COT  a 


Vu 


(A7) 


or 


co^  = CO  = - cos-1  (-a/y^ 


(A8) 


where  co^.  is  the  frequency  at  which  = -1  and  Z;  = 0,  for  given  values  of  t,  a,  and  A. 
Substituting  co^  from  Eq.  (A8)  for  co  in  Eqs.  (A6)  and  (A7),  we  obtain 

= cOp  =-^  cos-1  i^(xly/2J^.  (A9) 

Equation  (A9)  is  the  exact  stability  condition  from  which  can  be  found  the  value  of 
A = Ap  that  satisfies  Eq.  (A9)  for  given  values  of  a and  t.  The  system  will  then  be  stable 
if  the  parameters  are  chosen  such  that  A < A^.  However,  Eq.  (A9)  is  a transcedental  equa- 
tion that  is  difficult  to  solve  and  does  not  provide  convenient  interpretation.  Fortunately, 
for  the  range  of  system  parameters  in  which  we  are  interested,  2A,  which  allows  us 
to  neglect  a-  in  Eq.  (A9)  and  approximate  cos-1  it/2.  Then,  from  Eq.  (A9) 

we  have 


and 


Ap  ^ 7r2/2T2 


(AlO) 


CJ^  =4=  tt/t. 

Therefore,  for  the  system  to  be  stable,  the  system  parameters  should  be  chosen  such  that 

At2  < 7r2/2.  (All) 


The  numerical  results  presented  in  this  report  validate  the  accuracy  of  the  approxima- 
tions made  in  arriving  at  Eq.  (All). 
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